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The study of dynamical gluon mass generation at the level of Schwinger-Dyson equation involves 
a delicate interplay between various field-theoretic mechanisms The underlying local gauge in- 
variance remains intact by resorting to the well-known Schwinger mechanism, which is assumed 
to be realized by longitudinally coupled bound state poles, produced by the non-perturbative dy- 
namics of the theory. These poles are subsequently included into the Schwinger-Dyson equation 
of the gluon propagator through the three-gluon vertex, generating a non-vanishing gluon mass, 
which, however, is expressed in terms of divergent seagull integrals. In this talk we explain how 
such divergences can be eliminated completely by virtue of a characteristic identity, valid in di- 
mensional regularization. The ability to trigger this identity depends, in turn, on the details of the 
three-gluon vertex employed, and in particular, on the exact way the bound state poles are incor- 
porated. A concrete example of a vertex that triggers the aforementioned identity is constructed, 
the ensuing cancellation of all seagull divergences is explicitly demonstrated, and a finite gluon 
mass is obtained. Due to the multitude of conditions that must be simultaneously satisfied, this 
construction appears to be exclusively realized within the PT-BFM framework. The resulting 
system of integral equations gives rise to a gluon mass that displays power-law running and an 
effective charge which, due to the presence of the gluon mass, freezes in the infrared at a finite 
(non- vanishing) value. 
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1. Introduction 

The gluon is massless at the level of the fundamental QCD Lagrangian, and remains massless 
to all order in perturbation theory. However, as Cornwall argued in the early eighties [|T]], the 
non-perturbative QCD dynamics generate an effective, momentum-dependent mass for the gluons, 
without affecting the local SU (3)c invariance, which remains intact. 

Given that the gluon mass generation is a purely non-perturbative effect, the natural framework 
to study it, in the continuum, are the Schwinger-Dyson equations (SDEs) of the theory. At the level 
of the SDEs the generation of such a mass is associated with the existence of infrared finite solutions 
for the gluon propagator [|l|, ^. In covariant gauges, the gluon propagator, A^v(^), has the form 



(1.1) 



where B, denotes the gauge-fixing parameter, and ^^v{<j) = g^iv — '?^9v/^^- The scalar factor A(g^) 
is given by ^r^{q^) = + Kl{q^), where n^v(^) = P^iv('?)n(^^) is the gluon self-energy; the 
dimensionless vacuum polarization, T\.{q^), is defined as Yl{q^) = q^JI{q^). So, in general, one 
looks for solutions with A^' (0) > 0. Such solutions may be fitted by "massive" propagators of the 
form A^^{q^) = q^ + m^{q^); m^{q^) is not "hard", but depends non-trivially on the momentum 
transfer q^. 

The pinch technique (PT) propagator, usually denoted by A(^^) in the literature [|l], ^, is the 
ideal quantity to study in this context, because it is independent of the gauge-fixing parameter (i^). 
Therefore, any statement about its infrared behavior, and in particular the generation of a gluon 
mass, is bound to be free of gauge artefacts. In recent studies, however, the tendency has been to 
focus on the gluon propagator in a fixed gauge, such as the conventional Landau gauge {t, = 0), 
instead of the privileged Feynman gauge {E,q = 1) of the background field method (BFM) 
which, quite remarkably, reproduces automatically the PT results (third item in [Q]). The main 
reason for this less optimal choice is the need to compare meaningfully the SDE results with those 
obtained from lattice studies, which, almost exclusively, are carried out in the Landau gauge. 

What the latest large-volume lattice studies reveal is crystal clear: The gluon propagator of 
pure Yang-Mills is infrared finite, both in SU{2) and SU{2>) [^. Interestingly enough, these 
recent lattice findings, striking as they may be, do not constitute the first indication of this very 
characteristic behavior; several earlier simulations had found qualitatively similar results, even in 
gauges other than the Landau ( see, e.g., [Q]). 

The aforementioned lattice results, in addition to whatever modifications they may induce to 
other formalisms aspiring to describe the infrared sector of QCD, they present a serious challenge 
even to the practitioners of the gluon mass generation (albeit, a pleasant one). Indeed, the SDE 
analysis must be further refined, and freed of whatever major or minor theoretical shortcomings 
one has been accustomed to live with in the past. The purpose of this talk is to report recent 
progress in this direction, and in particular on the solution of the annoying problem of seagull 
divergences |^, which has afflicted this approach from the first day of its invention |[T]]. Turns out 
that the elimination of the seagull divergences brings about an additional advantage, namely the 
separation of the SDE into two coupled equations, furnishing uniquely the running of the effective 
charge and the gluon mass in the entire range of Euclidean momenta. 
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2. Seagull divergences: a perennial nuisance to gluon mass generation 

In order to obtain massive solutions gauge-invariantly, it is necessary to invoke the well-known 
Schwinger mechanism [|lO|]. The basic observation is that if, for some reason, Tl{q^) acquires a pole 
at zero momentum transfer, then the vector meson becomes massive, even if the gauge symmetry 
forbids a mass at the level of the fundamental Lagrangian. Indeed, it is clear that if the vacuum 
polarization Tl{q^) has a pole at = with positive residue /x^, i.e. = jJi^ /q^, then (in 

Euclidean space) (^^) = q^ + pL^. Thus, the vector meson becomes massive, A^' (0) = pL^, even 
though it is massless in the absence of interactions {g = 0). There is no physical principle which 
would preclude Yi{q^) from acquiring such a pole. In a strongly-coupled theory like QCD this may 
happen for purely dynamical reasons, since strong binding may generate zero-mass bound-state 



excitations [11|. The latter act like dynamical Nambu-Goldstone bosons, in the sense that they are 
massless, composite, and longitudinally coupled; but, at the same time, they differ from Nambu- 
Goldstone bosons as far as their origin is concerned: they do not originate from the spontaneous 
breaking of any global symmetry [[T]]. 

Of course, in order to obtain the full dynamics, such as, for example, the momentum-dependence 
of the dynamical mass, one must turn eventually to the SDE that governs the corresponding gauge- 
boson self-energy. The way the Schwinger mechanism is integrated into the SDE is through the 
form of the three-gluon vertex. The latter, even in the absence of mass generation, constitutes a 
central ingredient of the SDE, and plays a crucial role in enforcing the transversality of the gluon 
self-energy. Therefore, an important requirement for any self-consistent Ansatz used for that vertex 
is that it should satisfy the correct Ward identity (WI) of the PT-BFM formulation, namely 

q^f^^^=^-^'^{k + q)-^-^l{k). (2.1) 

In addition, in order to generate a dynamical mass one must assume that the vertex contains dy- 
namical poles, which will trigger the Schwinger mechanism when inserted into the SDE for the 
gluon self-energy. 

The point is that the full realization of the procedure outlined above is very subtle. In particular, 
even though the use of a three-gluon vertex containing massless poles and satisfying the correct WI 
leads indeed to a transverse and infrared finite self-energy (i.e. A^^(O) 7^ 0), as expected, the 
actual value of A^^(O) has always been expressed in terms of divergent integrals, of the form (see, 
e.g.,iii) 

A"i (0) = ci /" A(^) + C2I k^A^{k) , (2.2) 
Jk Jk 

where, in dimensional regularization (DR), = {2n)^^ Jd^k, with d = 4 — e the dimension 
of space-time. This is not a problem, in principle, provided that the divergent integrals appearing 
on the rhs of ( ^^ ) can be properly regulated and made finite, without introducing counterterms 
of the form ml{Aly)Aj^, which are forbidden by the local gauge invariance of the fundamental 
QCD Lagrangian. Various regularization procedures of increasing sophistication have been tried 
out over the years, but the resulting (regularized) A^' (0) remained theoretically ambiguous. Given 
how nicely all other pieces of the puzzle fit together, the underlying impression has always been 
that the "seagull problem" had to do with some (not fully understood) subtlety rather than some 
intrinsic "need" of the theory to produce quadratic divergences. 
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Figure 1: The "one-loop dressed" SDE for the photon self-energy. 



3. The "seagull identity": how to keep the photon massless (if you must) 

It is most instructive to understand what happens in a theory where the seagull terms do not ap- 
pear due to the self-interactions of a gauge boson that has acquired a mass dynamically, but rather 
because the theory has scalar particles that are massive at tree-level. These scalars interact with 
the gauge boson, and contribute seagull terms to its vacuum polarization. The question is: if the 
gauge boson must remain massless, how do the seagull contributions disappear from that vacuum 
polarization? To see how this happens, let us turn to scalar QED, where the aforementioned cir- 
cumstances (massive scalars, must-be massless photon) are realized, and study the SDE governing 
the photon self-energy. 

At the "one-loop dressed" level the SDE for the photon self-energy reads (Fig. [I]) 



n^v{q)=e'jTf&{k)&{k + q)Ty+e'jTfl&{k), 



(3.1) 



where ^^{k) is the fully-dressed propagator of the scalar field. Fy is the fully dressed photon- 
scalar vertex, whose tree-level expression is given by F^''^ = —i{lk + q)^. Moreover, the bare 
quadrihnear photon-scalar vertex is given by T^^l = lig^^. The photon-scalar vertex F^ and the 
scalar propagator ^ are related by the Abelian all-order WI 

q''ry = &-\k + q)-&-\k). (3.2) 

It is fairly easy to demonstrate that, by virtue of (3^), ^^n^v(^) = 0, and that Tl{q^) reads 

-2ie^ 



n{q') 



d-l 



&{k)&{k + q)k^T^ 



d J^S){k) 



(3.3) 



Let us compute from (33) the one-loop expression for n(^^), to be denoted by IT^^' ( 



n(^)(^^) 



-le 



d-l 



{Ak^-q')%{k)%{k + q)-2d / ^o(^) 
k Jk 



(3.4) 



where ^o(^) = [k^ — "^^) ^ ■ Taking the limit ^ — > 0, we find 



n(i)(o) 



-4/g^ 
d-\ 



%{k) 



(3.5) 



Of course, there is no doubt that the photon remains massless perturbatively, i.e. we must have 
that n(^)(0) = 0. However, the way this requirement is realized is rather subtle: the rhs of (R3) 
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vanishes indeed, by virtue of an identity that is exact in DR, namely 

k'^ d r I 



- m 



2^2 



d 



m 



2 ■ 



(3.6) 



Thus, the perturbative masslessness of the photon is explicitly realized and self-consistently en- 
forced within the DR. Eq.( |3.6| ) may be cast in a form that is particularly suggestive for the analysis 
that follows, namely 

,2d%ik) d 



(3.7) 



We now return to the general Eq.([3.3[). In order to analyze it further we must furnish some 
information about the form of . Of course, any meaningful Ansatz for must satisfy the WI of 
|), or else the transversality of n^v('?) will be compromised from the outset. The form obtained 



by Ball and Chiu [12], after "solving" the WI, under the additional requirement of not introducing 
kinematic singularities, is (we omit the identically conserved part of the vertex) 



{2k + q)^ 



[&-\k + q)-&-\k)\ . 



^ {k + qY-B 
This vertex, when substituted into ( |3.3D , yields 

"2 r r ^^{kJ^q)_^(k) 



(3.8) 



U{q' 



le 



d-l 



{Ae-q'Y- 
k [k + qY 



Taking the limit of Eq.([3.9|) as ^ — > 0, using that 

2{k + q)-2{k) 



{k + qY-k^ 



d^jk) 



-k^ 



+ 



+ 2d I &{k) 
k 



(3.9) 



we have that 



n(0) 



d-l 



,d&{k) 

'~dk2~ 



(3.10) 



(3.11) 



Of course, we must have that n(0) = 0, given that there is nothing in the dynamics that could 
possibly endow the photon with a mass; in particular, the Schwinger's mechanism is "switched 



off", i.e. we have not introduced dynamical poles, and, given the form of ( |3.8| ), neither kinematic 



ones, which might simulate the dynamical ones at the level of the SDE. Thus, the rhs of ( |3.1 1| ) must 
vanish, and therefore, we must have that 



dk'^ 



(3.12) 



which is the non-perturbative generalization of (3.7). 



Note a crucial point: the seagull terms appearing in ( [3.11| ) cannot be set to zero individually, 
because the scalar propagator inside them is massive (already at tree-level): the only way to keep 



the photon massless, is to employ (3.12), which cancels them against each other. For example, if 
the term Jj^ &{k) on the rhs were multiplied by any factor other than {d/2) one would be stuck with 
seagull divergences. 
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Figure 2: "One-loop dressed" gluonic graphs of the SDE for the PT-BFM gluon self-energy. 

4. Massive gluons: Schwinger mechanism and seagull identity in a delicate balance 

Turns out that the construction presented in the previous section generalizes in the context of 
pure Yang-Mills theory, such as quarkless QCD, but only within the PT-BFM formalism! As has 



been explained in detail in the recent literature [g, |13p, this latter formalism allows for a gauge- 
invariant truncation of the SD series, in the sense that it preserves manifestly and at every step 
the transversality of the gluon self-energy. Specifically, for the case at hand, we will consider 
only the "one-loop dressed" part of the gluon SDE that contains gluons, shown in Fig. ^, leaving 
out (gauge-invariantly!) the "one-loop dressed" ghost contributions and all "two-loop dressed" 
diagrams. Note that the Feynman rules used to build this SD series are those of the BFM [^]; in 
particular, the external gluons (distinguished by the grey circles attached to them) are treated as if 
they were background gluons. As we will see in a moment, the form of these vertices is crucial 
for obtaining from the SDE precisely the right combination of terms (and with the correct relative 
weights) that appears in ( [3.12| ). 



In order to reduce the algebraic complexity of the problem, we drop the longitudinal terms 
from the gluon propagators inside the integrals, i.e. we set A^^g — > —igf^pA. This does not compro- 
mise the transversality of n^v(^) provided that we do the same on the rhs of the WI satisfied by 
Fva/?! namely we have simply 

q^fy^p = [A-\k + q)-A-\k)]go,p , (4.1) 



instead of the full WI given in (^). Then, the SDE corresponding to Fig. g reduces to 



(4.2) 



where Ca the Casimir eigenvalue of the adjoint representation [Ca = A'^ for SU {N)\. The vertex 
T^^^^p[q^k,—k — q) is the bare three-gluon vertex in the Feynman gauge of the BFM, and T^f^^ 
denotes its fully-dressed version. 



The function A(^) appearing on the Ihs of ( |4.2[ ) is the scalar part of the gluon propagator in the 
BFM, i.e. two background gluons entering. A(^) is related to the standard A(^), defined in the 
gauges, by means of the powerful identity A(^)[l + G{q')Y = A{q), where G{q') is an auxiliary 



two-point function [14], which, quite remarkably, coincides in the Landau gauge with the well- 
known Kugo-Ojima function (see, e.g. talk of Daniele Binosi in these proceedings [^]). We will 
next set G{q') = 0, i.e. we effectively assume that, inside the quantum loops, A{q) = A{q). 
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Figure 3: The SDE for the three-gluon vertex. All kernels are one-particle irreducible, and the pole is 
not kinematic but dynamical (purely non-perturbative); physically it corresponds to a (composite) Goldstone 
mode, necessary for maintaining the local gauge invariance. 



At this point enters the new ingredient: a judicious Ansatz for the three-gluon vertex which, in 
addition to satisfying ( pTTl ) will allow us to use the seagull identity ( |3.12 ) and get a non- vanishing 
and finite A" ^(0). 

To begin with, let us first write A^\q) in the alternative form (in Minkowski space) 



(4.3) 



The tree-level result for A ^ (q) is recovered by setting H ^ (^) = 1 and = 0. Then, an appro- 
priate Ansatz for Ty^p is given by 

'{k + qfH-\k + q)- k^H- ^ (k) ' 



ir 



flap 



{k + qY-B 



(4.4) 



where the term V^^^g contains the non-perturbative contributions due to bound-state poles associ- 
ated with the Schwinger mechanism. Note that we must have 



q^y^iap = [m^{k)-m^{k + q)]gap , 



(4.5) 



in order for the r^„j} of Eq. ( [4.4| ) to satisfy (by construction) the correct WI of ( |4. 1[ ). 

The Ansatz of ( p] ) mimics that of Eq. ( |3.8| ) to the extent that the first term contains the 
right structure to produce, when inserted into the first term on the rhs of (12), the derivative term 
appearing on the Ihs of ( |3.12 ). The rhs of ( 3.12 ) is already there: it is the second term on the rhs of 
(pT^), originating directly from the seagull diagram (a2)- 

Similarly, a simple Ansatz for V^ap that captures the two essential characteristics of having a 
(composite), longitudinally coupled poles, and satisfying the WI of (4.5) is 



V, 



nap 



V a+V^ a 



where 



V, 



r,2 



n?{k)-m^{k + q) 



(4.6) 



(4.7) 



and with the transverse part Vj^^p satisfying q^V^^^p = 0. We can write the vertex of ( p^ equiva- 
lently as 

'A-^{k + q)-A-\k) 



{k + qY-B 



vap ^ fiaP' 



(4.8) 
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with 



where 



V 



nap — ^tiap + V^iap ' 



rrp-{k + q)-np-{k) 



?ap 



(4.9) 



(4.10) 



The term V^„^ is a residual piece, acting as an additional (non-perturbative) vertex term, originating 
from forcing the vertex to assume the form of (18). This last way of writing r^"/^ makes the use of 



the basic identity of Eq.( [3.12| ) immediate. Thus, after these rearrangements, we have that the final 
non-perturbative effective vertex V ^ccp must be transverse, q^V ^(j^p = 0. 

Substituting for the T^'^^ on the rhs the expression given in (^) we obtain after simple algebra 



A-\q' 



2{d-l] 



niq)+nff,{q) 



(4.11) 



with 



n(^) = {ld-S)q^ 



A{k + q)-A{k) 
\ {k + qf-k^ 



+ Ad 



A{k + q) 



{k + qf 



and 



n^(^) = fr^^ipA{k)A{k + q) [V' + {V' + r}rP , 



(4.12) 



(4.13) 



The term in square brackets on the rhs of (4.12) has exactly the structure needed for employing 
( |3.12 ). Note the perfect balance of relative coefficients required for this to happen! This becomes 
possible within the PT-BFM framework thanks to the special vertices appearing in the SDE; instead, 
in the conventional SD formulation (e.g., in the gauges) it would be very difficult to obtain the 
precise combination of terms needed for implementing ( 3.12| ). Evidently, by virtue of ( 3.12 ) it is 
clear that n(0) = 0. Thus, the part of the calculation determining n(^) is very similar to that of 
scalar QED, in the sense that it leads to total seagull annihilation, keeping the gluon massless. 

On the other hand, the term Ilff,{q), not present in the scalar QED study, makes it possible to 
have A^' (0) ^ for the gluons. Assuming, for simplicity, that the dominant contribution in ( |4.13 ) 
comes from V^, we obtain 



2d r 

^m{q) = T / tfiHk)A{k)A{k + q)[{k + qf-k^ 

q^ Jk 



(4.14) 



Now, in the limit — > (in Euclidean space) we have that 



lim i 4 / m^ik)A{k)A{k^ 



-q)[{k + qf 



-\ I k^A\k^)[fh\k^)'^ . (4.15) 
2 Jk, 



where the "prime" denotes differentiation with respect to k^. Note that a monotonically decreasing 
mass, [ffp-{k^)\' < 0, guarantees the positivity of m^(0) (in Euclidean space) 

An important consequence of this analysis is that Eq. ( [4.11| ) can be split unambiguously into 
two parts, one that vanishes as — and one that does not. In fact, using (1-.3) on the Ihs of 
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( |4.11[ ), we can assign the two types of contributions into two separate (but coupled) equations, 
namely (Minkowski space) 



g^H-\q)=q'-:^0^niq), (4.16) 

t^Hq) = 5^^n-(^)- (4.17) 

As we will see shortly, the first equation will determine the momentum dependence of the effective 
charge, and the second the running of the gluon mass. 

5. Effective charge and gluon mass: coupled but unique 

It is well-known that, due to the Abelian Wis of the PT-BFM Green's functions, the product 

forms a renormalization-group (RG)-invariant (/X -independent) quantity In order to reahze 



Eq.(5.1) non-perturbatively, first set 

m^{q^)=m\q^)H-\q^), (5.2) 
where m^{q^) is assumed to be the RG-invariant dynamical gluon mass. Then 

Aig')= "KLy (5-3) 
q^ +m^[q^) 

and from the requirement that g^A{q) must be RG-invariant we have that 

g^H{q^)=f{q^). (5.4) 
Therefore, we finally arrive at the RG-invariant combination 

diq^)^g^Aiq^)=fiq^)A{q^), (5.5) 

with 

Hq^)= 2^ \f 2v (5.6) 

q^ +m^[q^) 

So, d{q^) is written as the product of two RG-invariant quantities: the dimensionless running 
couphng g^{q^) and the dimensionful "massive" gluon propagator A{q^). 

We next cast our analysis in terms of the RG-invariant quantities defined above. The use of 
the spectral representation [|T], |T7|] for A{q^), namely 

A{q^)=jdX' (5.7) 
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Figure 4: Numerical solutions displaying the momentum-dependence of the effective charge and the gluon 
mass, for three different values of the renormalization point ji. 



results in a spectacular simplification, because it "solves" the combination {^j^^^jrS^} appearing 



in (4.12). After a series of standard assumptions one finally obtains 



1 



1 



and 



4z 



4^x1/2 



4z 



4z 



^^(1 + ^) (1-772) ^(^) 



1/2 



7.2(^2 



2b 

y 



^{q') / dyym\y)^{y) - 

2 Ja^ 



dyy^K\y)t{y)[m\y)^ 



(5.8) 



(5.9) 



where b = \QC\/A%k^; the discrepancy from the factor b=\ \Cp^/A^K^, namely the first coefficient 
of the QCD one-loop j3-function, is due to the (gauge-invariant!) omission of the ghost loops. 



To study the behavior of the solutions of (5^) for asymptotically large q^, set h{x) — > 1 /x and 
K{y) — > 1/j to arrive at 



m^{q^)\nq^ 



1 



dym^{y) 



1 



dyf{y)[mHy)]' 



(5.10) 



It is relatively straightforward to establish that the asymptotic solutions of ( 5.10D display power-law 
running. Indeed, substituting on both sides of (5.10) a of the form 



^2 



(5.11) 



one recognizes that the second term on the rhs of ( 5.11 ) is subleading, and that ( [5.1 1| ) is a solution 



of (|5.10|) provided that 7=1. 



We next solve numerically the two coupled integral equations, renormalizing them at three 
different points, namely /i = {4, 10,91}GeV, with a{ix^) = g^{n)/4n = {0.341,0.229,0.127}, 
respectively. In Fig. ^ we show the results for a{q^) and m^{q^); for either quantity we see clearly 
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RGI combination - d(q ) 

M = 91 GeVa(n^) = 0.127 

M = 1 GeV a(n") = 0.229 

H = 4 GeV a(n') = 0.341 




1E-3 0,01 



100 1000 10000 



Figure 5: The RG-invariant product d{q^) obtained by combining the results for Cic{q^) and m 



that the three curves merge practically into a single one, thus confirming numerically their jj.- 
independence, expected on formal grounds. The solutions for a{q^) may be fitted by the physically 
motivated functional form [p, namely 

1 



(5.12) 



with t = 3.7 and A = 645MeV. The behavior of m^{q^) in the entire range of momenta can be 
accurately described by the following parametrization 



tfln 



q'^+ml 



In 



A2 



In 



/(0,mi 



A2 



-3/5 



(5.13) 



where f{q^,ml) = p\m}. + p2 2'"° i » with pi = — 1/2, p2 = 5/2, and niQ = 612MeV. Finally, with 
the help of Eq. (5.5) we can construct d{q^) out of the numerical solutions for a{q^) and m^{q^); the 



result is shown in Fig. ^ Obviously, since d{q^) is built out of two quantities that are individually 
independent of jj., it too turns out to be /X -independent; this property is clearly observed in Fig. B. 



6. Conclusions 

The analysis presented here demonstrates that the appearance of seagull divergences in gluon 
mass generation is caused by a subtle mismatch between two parallel field theoretic mechanisms. 
Specifically, the Schwinger mechanism, which requires the appearance of massless poles in the 
three-gluon vertex, distorts the mechanism responsible for the cancellation of the seagull diver- 
gences, unless the poles enter into the gluon vertex in a very particular way. A concrete example 
of a vertex that does not produce any clash between these two mechanism has been given, and 
the implications for the resulting SDE have been worked out. The elimination of the seagull di- 
vergences allows the resulting SDE equation to be separated unambiguously into two distinct dy- 
namical equations, determining the gluon mass and the QCD effective charge. This, in turn, is 
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a significant improvement over the existing approaches, where the infrared behavior of these two 
quantities had to be extracted (not without a certain ambiguity) from d{q^). It is clear that the 
methodology outlined here should be applied to gauges where the results can be directly compared 
to lattice simulations (such as the Landau gauge). Given that the solutions for a{q^) and np-{q^), 
and therefore for d{q^) and/or A(^^), are expected to be unique, one should be able to test if the 
freezing value A(0) obtained within this new SDE treatment coincides with that seen on the lattice. 
We hope to report progress in this direction in the near future. 
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